Journal of Climate 

How well do global climate models simulate the variability of Atlantic tropical cyclones 

associated with ENSO? 



--Manuscript Draft- 

Manuscript Number: 


Full Title: 

How well do global climate models simulate the variability of Atlantic tropical cyclones 
associated with ENSO? 

Article Type: 

Article 

Corresponding Author: 

Hui Wang 

Climate Prediction Center 

Camp Springs, MD UNITED STATES 

Corresponding Author's Institution: 

Climate Prediction Center 

First Author: 

Hui Wang 

Order of Authors: 

Hui Wang 
Lindsey Long 
Arun Kumar 
Wanqiu Wang 
Jae-Kyung E. Schemm 
Ming Zhao 
Gabriel A. Vecchi 
Timorhy E. LaRow 
Young-Kwon Lim 
Siegfried D. Schubert 
Daniel A. Shaevitz 
Suzana J. Camargo 
Naomi Henderson 
Daehyun Kim 
Jeffrey A. Jonas 
Kevin J. E. Walsh 

Abstract: 

The variability of Atlantic tropical cyclones (TCs) associated with El Nino-Southern 
Oscillation (ENSO) in model simulations is assessed and compared with observations. 
The model experiments are 28-yr simulations forced with the observed sea surface 
temperature from 1982 to 2009. The simulations were coordinated by the U.S. 

CLIVAR Hurricane Working Group and conducted with five global climate models 
(GCMs) with a total of 16 ensemble members. The model performance is evaluated 
based on both individual model ensemble means and multi-model ensemble mean. 

The latter has the highest anomaly correlation (0.86) for the interannual variability of 
TCs. Previous observational studies show a strong association between ENSO and 
Atlantic TC activity, as well as distinctions in the TC activities during eastern Pacific 
(EP) and central Pacific (CP) El Nino events. The analysis of track density and TC 
origin indicates that each model has different mean biases. Overall, the GCMs 
simulate the variability of Atlantic TCs well with weaker activity during EP El Nino and 
stronger activity during La Nina. For CP El Nino, there is a slight increase in the 
number of TCs as compared with EP El Nino. However, the spatial distribution of track 
density and TC origin is less consistent among the models. Particularly, there is no 
indication of increasing TC activity over the U.S. southeast coastal region as in 
observations. The difference between the models and observations is likely due to the 
bias of vertical wind shear in response to the shift of tropical heating associated with 
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ABSTRACT 


The variability of Atlantic tropical cyclones (TCs) associated with El Nino-Southern 
Oscillation (ENSO) in model simulations is assessed and compared with observations. The 
model experiments are 28-yr simulations forced with the observed sea surface temperature from 
1982 to 2009. The simulations were coordinated by the U.S. CLIVAR Hurricane Working 
Group and conducted with five global climate models (GCMs) with a total of 16 ensemble 
members. The model performance is evaluated based on both individual model ensemble means 
and multi-model ensemble mean. The latter has the highest anomaly correlation (0.86) for the 
interannual variability of TCs. Previous observational studies show a strong association between 
ENSO and Atlantic TC activity, as well as distinctions in the TC activities during eastern Pacific 
(EP) and central Pacific (CP) El Nino events. The analysis of track density and TC origin 
indicates that each model has different mean biases. Overall, the GCMs simulate the variability 
of Atlantic TCs well with weaker activity during EP El Nino and stronger activity during La 
Nina. For CP El Nino, there is a slight increase in the number of TCs as compared with EP El 
Nino. However, the spatial distribution of track density and TC origin is less consistent among 
the models. Particularly, there is no indication of increasing TC activity over the U.S. southeast 
coastal region as in observations. The difference between the models and observations is likely 
due to the bias of vertical wind shear in response to the shift of tropical heating associated with 
CP El Nino, as well as the model bias in the mean circulation. 
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1. Introduction 


It is well known that El Nino-Southern Oscillation (ENSO) strongly influences the 
interannual variability of Atlantic tropical cyclones (TCs). El Nino (La Nina) tends to suppress 
(enhance) Atlantic seasonal TC activity (e.g., Gray 1984; Pielke and Landsea 1999; Bell and 
Chelliah 2006). Although other climate modes, such as the Atlantic Meridional Mode, the North 
Atlantic Oscillation, and the Madden Julian Oscillation, also modulate North Atlantic tropical 
cyclone activity (e.g. Kossin et al. 2010), here our focus will be placed on ENSO. The state of 
ENSO is one of the key climate factors considered by the National Oceanic and Atmospheric 
Administration (NOAA) for their Atlantic hurricane season outlooks (NOAA 2013). 

Using observational data, Kim et al. (2009) found distinct differences in Atlantic TC 
activity associated with eastern Pacific (EP) El Nino and central Pacific (CP) El Nino. EP El 
Nino is the conventional El Nino with the warmest sea surface temperature (SST) anomalies in 
the tropical eastern Pacific, whereas CP El Nino or El Nino Modoki (Ashok et al. 2007) is a non- 
conventional El Nino with the warmest SST anomalies in the tropical central Pacific. The zonal 
shift of the warm SST anomalies indicates a change in tropical heating and consequent changes 
in atmospheric response. 

A composite analysis of TC track density anomaly in Kim et al. (2009, their Fig. 2) 
displays coherent weakening in TC activity over the Caribbean Sea, Gulf of Mexico, and U.S. 
Atlantic east coast region during EP El Nino and strengthened TC activity over the same regions 
during La Nina. Surprisingly, the composite for CP El Nino is also opposite to that for EP El 
Nino over these regions and closely resembles the La Nina composite. The results suggest a 
higher chance of landfalling TCs along the Gulf coast and U.S. east coast during CP El Nino 
than during EP El Nino. 
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It is well recognized that global climate models (GCMs), even at a low resolution, are 
able to simulate the interannual response of North Atlantic TCs to ENSO (e.g. Camargo et al. 
2005, Zhao et al. 2009). However, given the distinctions in the Atlantic TC activity associated 
with different El Nino types revealed in observations (Kim et al. 2009), it is interesting to know 
whether state-of-the-art GCMs can reproduce the different response to the two types of El Nino. 
Such a model capability in distinguishing the responses of Atlantic TCs to different ENSO 
patterns is also important to both dynamical (e.g., Schemm and Long 2009) and statistical- 
dynamical (e.g., Wang et al. 2009; Vecchi et al. 2011) hurricane seasonal prediction systems. 

With a primary focus on climate modeling studies of TCs, the U.S. Climate Variability 
and Predictability Research Program (CLIVAR) launched a Hurricane Working Group (HWG) 
in 2011 (U.S. CLIVAR 2011). To improve understanding of the interannual variability and 
trends in TC activity, as well as projections of future TC activity under a warming climate, the 
HWG initiated a series of simulations with high-resolution atmospheric GCMs (Walsh et al. 
2013). One set of simulations is the interannual experiment which is an Atmospheric Model 
Intercomparison Project (AMIP) type of simulations with multiple GCMs and forced with the 
same observed time-varying SST from 1982 to 2009. This set of simulations provides necessary 
data to characterize TC response to ENSO in climate models. 

This study aims to evaluate the performance of high-resolution GCMs in simulating the 
interannual variability of Atlantic TCs associated with ENSO. The assessment is based on the 
analysis of AMIP-type simulations with five GCMs and comparisons with observations. The 
analysis is in collaboration with HWG members to target one of the HWG objectives involving 
improved understanding of interannual variability of TC activity. The following three scientific 
questions are to be addressed in this study. How is the overall performance of GCMs in 
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simulating the variability of Atlantic TCs? What are the characteristics of Atlantic TCs 
associated with ENSO in the models? What are the possible explanations for the differences 
between the models and observations? The study is expected to provide some insights into the 
basic characteristics of Atlantic TC activity associated with different types of ENSO in GCMs. 

This paper is organized as follows. Section 2 provides a brief description of data, 
models, and analysis methods used. Section 3 characterizes the Atlantic TC activity associated 
with ENSO in observations. The performance of GCMs in simulating the variability of the 
Atlantic TCs is assessed in section 4. Some possible explanations for the differences between 
the model simulations and observations are explored in section 5. Conclusions are given in 
section 6. 

2. Data and models 

The data used in this study consist of SST, Atlantic TC tracks, precipitation, 200-hPa and 
850-hPa zonal winds over a 28-yr (1982-2009) period from both observations and simulations 
with five atmospheric GCMs. For observations, the SST data are taken from the Hadley Centre 
Sea Ice and Sea Surface Temperature (HadlSST) data set (Rayner et al. 2003) on a 1° x 1° 
(latitude x longitude) grid. The 28-yr monthly mean SSTs were also prescribed as low boundary 
forcing for the GCMs. The Atlantic TC track data are from the National Hurricane Center 
Atlantic Hurricane Best Track Data (HURDAT2; Landsea et al. 2004). The precipitation data 
are from the Climate Prediction Center (CPC) Merged Analysis of Precipitation (CMAP) data set 
(Xie and Arkin 1997). The 200-hPa and 800-hPa zonal winds used to derive vertical wind shear 
are from the National Centers for Environmental Prediction - Department of Energy (NCEP - 
DOE) reanalysis 2 (R2; Kanamitsu et al. 2002). Both precipitation and zonal winds are monthly 
mean data on a 2.5° x 2.5° grid. 
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The five GCMs employed for the HWG interannual experiments (1982-2009) are the 
Florida State University (FSU) model (Cocke and LaRow 2000), Geophysical Fluid Dynamics 
Laboratory (GFDL) model (Zhao et al. 2009), National Aeronautics and Space Administration 
(NASA) Goddard Institute for Space Studies (GISS) model E2 (Schmidt et al. 2013), NASA 
Goddard Space Flight Center (GSFC) GEOS-5 model (Rienecker et al. 2008; Molod et al. 2012), 
and NCEP Global Forecast System (GFS) model (Saha et al. 2013). Table 1 lists the number of 
ensemble runs, model data resolutions, which are close to model resolutions, and TC tracking 
algorithms for the five models. The ensemble members vary from two to five with a total of 16 
runs. Horizontal resolutions range from about 0.5° to 1°. TC track data were provided by each 
modeling group with different tracking algorithms. More detailed descriptions of the models can 
be found in Walsh et al. (2013). 

The Atlantic TC activity is quantified by the annual total number of TCs, as well as the 
spatial distribution of track density and TC origin. Given the spatially discrete nature of TC 
tracks, the track density is derived as follows: (a) the number of TCs passing through each 5° x 
5° box analyzed on a 1° x 1° grid resolution during an entire hurricane season is first counted; 
and (b) the TC counts are then averaged with the TC numbers in the 5° x 5° boxes for eight 
surrounding grid points with a weighting coefficient of 0.5 for the center grid point and 1/16 for 
each surrounding grid point. This is done in the same way as Kim et al. (2009) to ensure a 
spatially smoothed distribution. Composites of SST, precipitation, and vertical wind shear 
anomalies averaged over August-October (ASO), the peak of the Atlantic hurricane season, are 
examined for different ENSO categories. The statistical significance of the composite anomalies 
is estimated by the Monte Carlo technique (e.g., Wilks 1995). The analysis is performed for both 
observations and multi-model ensemble (MME) mean, as well as individual model ensemble 
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means. The MME mean is obtained by averaging individual model ensemble means. In this 
way, each model is treated with an equal weight for the MME, regardless of the number of 
ensemble members. 

3. Variability of Atlantic TCs associated with ENSO in observations 

During the 28-yr period (1982-2009), there were five EP El Nino (1982, 1986, 1991, 

1997, and 2006) and five CP El Nino (1987, 1994, 2002, 2004, and 2009) years identified based 
on the definition of McPhaden et al. (2011), and eight La Nina years (1983, 1984, 1988, 1995, 

1998, 1999, 2005, 2007). Figure 1 shows the composite of ASO seasonal mean SST anomalies 
for EP El Nino, CP El Nino, and La Nina, respectively. Compared to EP El Nino (Fig. la), the 
SST anomalies in CP El Nino (Fig. lb) shift towards the west. This may lead to significant 
changes in tropical heating for the atmosphere between the two types of El Nino. The amplitude 
of the CP El Nino SST anomalies (~ 1 K) is smaller than the EP El Nino (~ 1.5 K), but 
comparable to the La Nina (~ 1 K, Fig. lc). 

Similar composites are shown in Fig. 2 for TC track density (top row) and track density 
anomaly (middle row), respectively, associated with the three ENSO types. In La Nina years 
(Fig. 2c), track density displays high values (> 1) across the North Atlantic basin. Areas with 
track densities greater than 1.5 are found in the central main development region (MDR; 10° - 
20°N, 20° - 80°W), the Gulf of Mexico, and U.S. east coastal region. In contrast, track density is 
relatively low over these regions for EP El Nino (Fig. 2a), but increases considerably for CP El 
Nino (Fig. 2b), particularly in the MDR and U.S. southeast coastal region. 

Consistent with the track density patterns, track density anomalies are generally below 
normal across the basin for EP El Nino (Fig. 2d), with the largest negative anomalies over the 
Gulf and MDR, and above normal during La Nina (Fig. 2f). Associated with CP El Nino (Fig. 
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2 e), positive track density anomalies are found over the MDR, the Caribbean Sea, Gulf coast and 
the southeast coast, and negative anomalies further to the east, as well as in the west Gulf of 
Mexico. The results indicate that relative to EP El Nino, there is a high chance of landfalling 
TCs along the U.S. southeast coast during CP El Nino. 

The spatial distributions of total TC origins for the three ENSO categories are also shown 
in Fig. 2 (bottom row). For a fair comparison with five EP El Nino and five CP El Nino, TC 
origins for La Nina are also shown for five episodes that occurred in the most recent years. 
There are increased TC origins over the MDR during CP El Nino (Fig. 2h) as compared to EP El 
Nino (Fig. 2g) and an additional increase of TC formation over the Gulf of Mexico during La 
Nina (Fig. 2i). 

Although the sample size for ENSO composites is very limited over the 28 years, the 
composite anomalies in Fig. 2 (middle row) are statistically significant above the 90% level. The 
anomaly patterns also resemble those in Kim et al. (2009) with longer records (57 yrs, 1950— 
2006). Additionally, the sampling issue can be partially addressed by using HWG interannual 
experiments which provide more atmospheric realizations than for the observations. Although 
the AMIP type of simulations does not increase the sample size of ENSO events, the ensemble of 
AMIP runs presented in the next section increases the sample size of atmospheric realizations for 
a fixed set of ENSO events. This can effectively enhance the signal-to-noise ratio (Kumar and 
Hoerling 1995) and thereby provide a more reliable estimate for the ENSO-forced variability of 
the Atlantic TCs. 

4. Variability of Atlantic TCs associated with ENSO in GCMs 

The climatology and interannual variability of the annual number of Atlantic TCs are 
examined first. Figure 3 shows the time series of the annual number of Atlantic TCs from 1982 
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to 2009 for both observations and model simulations, including MME mean and individual 
model ensemble means. Both observations and MME display an upward trend over the 28 -yr 
period. The grey shading in Fig. 3 denotes the range of ± one standard deviation of the spreads 
of the five individual model ensemble means around the MME mean. Over 80% (23 out of 28 
yrs) of the observations fall into this range. Obviously, the GFS model has very high numbers of 
TCs and the GISS model has low numbers of TCs. 

Table 2 summarizes the TC statistics for the observations and model simulations, 
including the climatological mean value, variance of interannual variability, linear trend over the 
28 years, anomaly correlation (AC) between the models and observations, and root-mean square 
error (RMSE). The GFDL model (12.7) and GSFC model (10.9) have a mean value close to the 
observations (1 1.7). In contrast, the climatology in the GISS model (6.2) is only about a half of 
the observations while the GFS model (22.0) has double the number in observations. The 
strength of the interannual variability in the GSFC and GFS models is comparable to 
observations and weaker in the other models and the MME. The linear trends in all models (~ 2 
TCs per decade) are weaker than in the observations (~ 4 TCs per decade). AC is highest for the 
MME (0.86), followed by the GFDL (0.74) and GFS (0.73) models. This implies that 74% of 
the observed interannual TC variance is captured by the time series of the MME mean number of 
TCs and 54% is captured by the GFDL and GFS models. Additionally, the MME has the 
smallest RMSE. Due to the large mean biases, the GFS and GISS models have relatively large 
RMSEs. In terms of the five parameters in Table 2 (i.e., mean, interannual variability, trend, AC, 
and RMSE), the overall performance of the MME, GFDL and GSFC models is better than that of 
the FSU, GISS, and GFS models. It should be noted that both the GFDL and GSFC models have 


8 


230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 


a higher resolution than the other three models. This may suggest that a GCM with a higher 
resolution gets better performance in simulating the interannual variability of Atlantic TCs. 

The average number of TCs for each ENSO category is examined in Table 3 and 
compared with the corresponding 28-yr climatology for both observations and simulations. In 
the observations, there are about 7, 10, and 15 TCs each hurricane season in EP El Nino, CP El 
Nino, and La Nina, respectively, equivalent to 58%, 87%, and 125% of the mean value (11.7). 
All models show consistent increases in the number of TCs from EP El Nino to CP El Nino and 
further increases to La Nina, except for the GSFC model. However, the changes in TC counts 
from one ENSO type to another in the models are much more conservative than in the 
observations. In the MME, for instance, there is a 15% increase in TCs from EP El Nino to CP 
El Nino and an additional 16% increase to La Nina in terms of the mean value. The 
corresponding changes in observations are 29% and 38%. The results indicate a weaker 
interannual variability of Atlantic TCs in the model simulations. It should also be noted that the 
MME mean approach may reduce the variability of TC counts in the models. 

The spatial characteristics of mean TC activity are presented in Fig. 4 for both 
observations and simulations in the form of 28-yr mean track densities and total TC origins 
during the entire 28 years. Compared to the observations (Fig. 4a), each model has different 
mean biases. Among the five models, the GFDL model (Fig. 4d) is closest to the observations 
for both the magnitude and spatial coverage of track density. The FSU, GSFC, and GFS models 
(Figs. 4c, 4f, 4g) have a very high track density (> 3) over the west MDR, east-central MDR, 
and most of the North Atlantic basin, respectively, whereas the GISS model (Fig. 4e) has a very 
low track density over the basin. The MME mean pattern (Fig. 4h) shows a higher track density 
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in the MDR than the observations (Fig. 4a) and a lower track density over the U.S. east coastal 
regions. Overall, the MME is better than most individual models. 

The TC origins in observations (Fig. 4b) are characterized by two regions with large 
populations, one over the MDR and the other over the Gulf of Mexico and adjacent sectors of the 
Atlantic Ocean and Caribbean Sea. The FSU, GSFC, and GFS models exhibit very dense TC 
origins over the central and to the south of the MDR (Fig. 4i), to the south of the east MDR (Fig. 
41), and to the south and east of the MDR (Fig. 4m), respectively. The GISS model shows a lack 
of TC formations over the east MDR. The GFDL model (Fig. 4j) and MME (Fig. 4n) have a 
distribution of TC origins closer to the observations than the other models. The model biases in 
the distribution of TC origins are consistent with the biases of track density and mean number of 
TCs. For example, the dense TC origins in the FSU and GSFC models (Figs. 4i and 41) lead to 
high track density over the regions to the northwest of the TC origins (Figs. 4c and 4f). If the 
unrealistic TC origins to the east of the MDR in the GFS model (Fig. 4m) are removed, the mean 
number of TCs is significantly reduced from 22.0 to 11.7, matching the observed value, and 
leading to a track density distribution much closer to the observations (not shown). 

Similar to the ENSO composites of track density for observations (Fig. 2, top row), Fig. 5 
displays the ENSO composites of track density for individual model ensemble means, as well as 
MME mean. In spite of the distinct biases in each model revealed in Fig. 4, the composites 
consistently show relatively low track densities during EP El Nino (left column) in all models 
and high track densities during La Nina in most models (right column), except for the GSFC 
model. Furthermore, there is a clear increase in track density from EP El Nino to CP El Nino 
(Fig. 5, middle column). 
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The corresponding composites for track density anomaly are illustrated in Fig. 6. The 
track density anomalies in the GCMs are generally below normal across the basin during EP El 
Nino (left column) and above normal during La Nina (right column). In some spots, the negative 
anomalies associated with EP El Nino (left column) become positive during CP El Nino (middle 
column). The results in Figs. 5 and 6 suggest that the GCMs are able to capture some of the 
observed features of the Atlantic TC activity associated with ENSO. Qualitatively, there is less 
TC activity associated with EP El Nino, more activity associated with La Nina, and increasing 
TC activity during CP El Nino with respect to EP El Nino. However, the patterns of track 
density vary from model to model and differ from observations. Particularly, there are no 
indications of increasing landfalling TCs along the U.S. southeast coast during CP El Nino in the 
model simulations. 

The modeled TC origins over five years from one ensemble member of each model are 
shown in Fig. 7 for each ENSO category. Relative to EP El Nino (left column), there are 
increases in the formation of TCs over or near the MDR during CP El Nino (middle column) and 
La Nina (right column) in some models, such as the GSFC and GFS models. Only the GFDL 
model shows some increase in TC origins at high latitudes between 20°N and 40°N, especially 
during CP El Nino. Unlike observations (Fig. 2i), there are no increases in TC origins over the 
Gulf of Mexico and west Caribbean Sea in all models during La Nina. This may be related to the 
model bias in simulating the TC formations over these regions (Fig. 4). The differences in TC 
origins among the three ENSO categories in the MME (Fig. 7, bottom row) are not as large as in 
the observations (Fig. 2, bottom row). This is another indication of relatively weak interannual 
variability of Atlantic TCs in GCMs. 


11 


297 

298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 


5. Possible explanations for model biases 

The changes in both the mean and variability of Atlantic TCs is accompanied by changes 
in atmospheric circulation (e.g., Goldenberg and Shapiro 1996; Goldenberg et al. 2001). 
Therefore, in order to understand the mean biases of TC activity in GCMs, Fig. 8 shows the ASO 
season climatology of vertical shear of zonal wind between 200 and 850 hPa derived from 
observations and mean biases for individual model ensemble means and the MME mean. The 
regions of weak mean vertical wind shear (< 10 m s' 1 , Fig. 8a) coincide with the regions of high 
mean track density and TC origins in observations (Figs. 4a and 4b). 

The mean bias in the vertical wind shear may account for the mean bias in Atlantic TC 
activity in some models. In the FSU model (Fig. 8b), for instance, a large negative bias of 
vertical wind shear (over -10 m s' 1 ) in the west MDR leads to a close-to-zero mean state of 
vertical wind shear, which favors the generation and development of TCs. This is consistent 
with the mean bias of high track density and TC origins over this region (Figs. 4c and 4i). In the 
GISS model (Fig. 8d), a positive bias of vertical wind shear in the east MDR enhances the mean 
vertical wind shear and prevents TCs from occurring over this area. As a result, TC tracks and 
TC origins shift towards the west (Figs. 4e and 4k). 

Both individual model ensemble means (Figs. 8b— 8f) and the MME mean (Fig. 8g) 
exhibit negative biases in vertical wind shear over and/or near the MDR and positive biases to 
the north, especially over the Gulf coast and U.S. southeast coast. Consequently, there are biases 
of high track density and dense TC origins at low latitudes and low track density and sparse TC 
origins over the Gulf and U.S. southeast coast in the models (Fig. 4). 

Figure 9 displays the composites of ASO season vertical wind shear anomalies associated 
with the three ENSO categories for observations (top row) and MME (bottom row), respectively. 
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Overall, the model circulation response to different ENSO SST anomalies agree with the 
observations, both with positive vertical wind shear anomalies to the south of 20 °N associated 
with EP El Nino (left column) and negative anomalies associated with La Nina (right column). 
The circulation response to CP El Nino is less significant or spatially coherent over the 
subtropical North Atlantic (middle column). This is likely due to the weak amplitude and small 
area-coverage of the CP El Nino SST anomalies (Fig. 1). Thus the atmospheric response may be 
weak (e.g., Wang et al. 2013). In spite of that, it is still evident that wind shear anomalies over 
the MDR are largely reduced as compared to EP El Nino, a condition that is more favorable for 
TC activity during CP El Nino. The results present in Fig. 9 are also consistent with the better 
simulations of Atlantic TC activity in GCMs for EP El Nino and La Nina than for CP El Nino. 

ENSO influences the Atlantic TC activity by altering vertical wind shear over the MDR 
through atmospheric teleconnection (e.g., Goldenberg and Shapiro, 1996). It may also change 
tropical Atlantic SST via local air-sea interaction (Enfield and Mayer 1997), which in turn 
affects the TC activity (Goldenberg et al. 2001). The composites of SST anomalies in Fig. 1 
suggest very weak Atlantic SST anomalies associated with ENSO in ASO. Furthermore, 
diagnostics of the ENSO modulation of TC activity using a genesis potential index identified 
vertical wind shear as one of the main environmental factors responsible for this modulation in 
the North Atlantic (Camargo et al. 2007). Therefore, the atmospheric response to tropical 
heating related to ENSO SST anomalies and atmospheric teleconnection are likely the primary 
processes responsible for the ENSO impact. 

The westward shift of warm SST anomalies from EP El Nino to CP El Nino (Fig. 1) may 
lead to changes in tropical heating. In the tropics, precipitation associated with deep convection 
is a good indicator of tropical heating in the atmosphere. Similar to Wang et al. (2012), the 
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composites of ASO season precipitation anomalies over the tropical Pacific are used to illustrate 
and verify the changes in tropical heating, as shown in Fig. 10. In both observations and the 
MME mean of the GCM simulations, associated with EP El Nino (Figs. 10a and lOd), there are 
positive precipitation anomalies across the central and eastern equatorial Pacific. Associated 
with CP El Nino (Figs. 10b and lOe), precipitation anomalies shift towards the west with no 
large anomalies over the eastern Pacific. In La Nina, negative precipitation anomalies cross the 
tropical Pacific (Figs. 10c and lOf). In general, the GCMs reproduce the observed major features 
of precipitation anomalies over the tropical Pacific for different types of ENSO. On the other 
hand, precipitation response to ENSO over the tropical North Atlantic (not shown) varies 
considerably across the models, which may contribute to the model diversity in simulating the 
TC variability associated with ENSO. 

There are also differences in precipitation between observations and simulations over the 
tropical Pacific, such as weaker precipitation anomalies in the models between 160°E and the 
dateline for all ENSO categories. These differences may be related to model convection 
schemes and model sensitivity to SST. Together with model biases in mean circulation (not 
shown), they may modify the Rossby-wave source (Sardeshmukh and Hoskins 1998) and thus 
affect the detailed structure of circulation response to ENSO. 

Figure 1 1 gives a simple example of changes in vertical wind shear associated with a 
westward shift of warm SST anomalies from the Nino-3 region (5°S - 5°N, 90° - 150°W) to the 
Nino-4 region (5°S - 5°N, 150°W - 160°E). First, the ASO season vertical wind shear anomalies 
are regressed against the Nino-4 and Nino-3 SST indices, separately. The differences between 
the two sets of regression coefficients are shown for observations (left panel) and the MME 
(right panel), respectively. Both the observations and the MME exhibit a similar large-scale 
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wave train pattern originating from the western and central equatorial Pacific and along a great 
circle route to tropical Atlantic. A close inspection of Fig. 1 1 reveals some differences in the 
changes of vertical wind shear over the tropical North Atlantic between the observations and 
simulations. Negative wind shear anomalies are found to the north of the MDR in the 
observations (Fig. 11a) whereas positive anomalies are found over the MDR in the MME (Fig. 
1 lb). The results illustrate the difference between the observations and GCMs in North Atlantic 
vertical wind shear response to the shift of tropical Pacific SST anomalies. The difference may 
cause further changes in the responses of Atlantic TCs to the shift of SST anomalies. 

6. Summary and conclusions 

Based on the analysis of the FIWG interannual experiments, the GCM’s performance in 
simulating the variability of Atlantic TCs associated with ENSO are assessed. The results 
indicate that each model has different mean biases in terms of track density and TC origin. 
Among the five models, the GFDL model with a relatively high resolution has the best 
performance. The MME mean has the highest anomaly correlation for the number of TCs and 
the least RMSE. Therefore, using an MME should be considered a better approach for 
dynamical hurricane season prediction than using a single model. Overall, the GCMs simulate 
the variability of Atlantic TCs well with weaker activities during EP El Nino and stronger 
activities during La Nina. For CP El Nino, there is a slight increase in the number of TCs as 
compared with EP El Nino. However, the spatial distribution of track density and TC origin is 
less consistent among the models. Particularly, there is no indication of increasing TC activity 
over the U.S. southeast coastal region as found in observations. The differences between the 
models and the observations may be due to the bias of vertical wind shear in response to the shift 
of tropical heating associated with CP El Nino, as well as the model bias in the mean circulation. 
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It should also be noted that there are limited sample sizes for both EP and CP El Nino events in 
the observations. The differences between EP and CP El Nino composites may not be just due to 
ENSO response, but also contain some random component. 

There are at least two factors that may affect the results presented in this paper. One is 
the model sensitivity to different SST data sets (e.g., LaRow 2013). For example, the FSU 
model forced with the NOAA optimum interpolation SST version 2 (OISST v2; Reynolds et al. 
2002) may improve the simulations of Atlantic TC activity with a better TC climatology (1 1.5) 
and RMSE (4.5) than those forced with HadlSST (Table 2). Knowledge of the model sensitivity 
to SST forcing may help estimate the uncertainty of the model simulated TCs. In this study, 
different TC tracking algorithms were employed by the five modeling groups for their GCMs 
(Table 1). Track density and TC origin in the models may also be sensitive to the algorithms 
used (e.g., Horn et al. 2013). A unified tracking algorithm may be helpful to reduce the related 
uncertainty for model assessment. 

The impact of ENSO on Atlantic TC activity may have some implications for projections 
of future TC variability under a warming climate. Studies have shown an increase in tropical 
Atlantic wind shear (Vecchi and Soden, 2007) and a reduction of Atlantic TCs associated with 
global warming with a high-resolution GCM (Zhao and Held 2010). In more recent studies, no 
robust changes in North Atlantic TC activity were found in the 21 st century simulations with 
low-resolution models (Camargo 2013; Tory et al. 2013). On the other hand, downscaling 
studies of these simulations lead to contradictory results, varying from a significant decrease 
(Knutson et al. 2013), ambiguous trends (Villarini and Vecchi 2013), to a significant increase 
(Emanuel 2013) in North Atlantic TC activity by the end of the 21 st century. In addition to 
possible changes in the mean TC activity, the variability of TC activity is also expected to 
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change as the intensity of CP El Nino (EP El Nino) would increase (decrease) under a warming 
climate (Kim and Yu 2012). In fact, CP El Nino has been documented to occur more frequently 
in the most recent two decades (Yeh et al. 2009), which could be a manifestation of global 
warming in observations. 

There is a possibility that the relationship between Atlantic TC activity and ENSO under 
the present-day climate found in Kim et al. (2009) might not be maintained under a warming 
climate. Indeed, changes in atmospheric teleconnection in response to ENSO have been detected 
in model simulations for the 21 st century (e.g., Stevenson 2012). This would add additional 
uncertainty to the future projection of Atlantic TC variability. Nevertheless, this study indicates 
the feasibility of utilizing high-resolution GCMs to assess the Atlantic TC activity associated 
with ENSO for climate change projections. 
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Figure Captions 

Fig. 1. Composites of ASO seasonal mean SST anomalies (unit: K) for (a) EP El Nino 
(1982, 1986, 1991, 1997, and 2006), (b) CP El Nino (1987, 1994, 2002, 2004, and 2009), and (c) 
La Nina (1983, 1984, 1988, 1995, 1998, 1999, 2005, and 2007) during 1982-2009. The 
anomalies circled by light lines are above the 99% significance level estimated by the Monte 
Carlo tests. 

Fig. 2. Composites of TC track density (top row) and track density anomaly (middle row) 
for (a),(d) EP El Nino, (b),(e) CP El Nino, and (c),(f) La Nina years, and distribution of TC 
origins during (g) five EP El Nino, (h) five CP El Nino, and (i) five La Nina years derived from 
observations. The anomalies circled by light lines (middle row) are above the 90% significance 
level estimated by the Monte Carlo tests. The boxes with dash lines denote the main 
development region (MDR; 10°-20°N, 20°-80°W). 

Fig. 3. Time series of annual number of Atlantic TCs from 1982 to 2009 for observations 
(OBS) and multi-model ensemble (MME) mean (thick lines with open circles), as well as 
individual model ensemble means (thin lines). Grey shading denotes the range of± one standard 
deviation of the spreads of the five individual model ensemble means around the MME mean. 

Fig. 4. Climatology of track density for (a) observations, (c)-(g) individual model 
ensemble means, and (h) MME mean, and 28-yr total TC origins for (b) observations, (i)-(m) 
one ensemble member of each model, and (n) MME total from one member of each model. The 
boxes with dash lines denote the MDR. 

Fig. 5. Composites of track density during EP El Nino (left column), CP El Nino (middle 
column), and La Nina (right column) for five individual model ensemble mean (top five rows) 
and MME mean (bottom row). The boxes with dash lines denote the MDR. 
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Fig. 6. Composites of track density anomaly during EP El Nino (left column), CP El Nino 
(middle column), and La Nina (right column) for five individual model ensemble mean (top five 
rows) and MME mean (bottom row). The anomalies circled by light lines are above the 90% 
significance level estimated. The boxes with dash lines denote the MDR. 

Fig. 7. Distribution of TC origins during five EP El Nino (left column), five CP El Nino 
(middle column), and five La Nina (right column) years from one ensemble member of each 
model (top five rows) and MME total from one member of each model (bottom row). The boxes 
with dash lines denote the MDR. 

Fig. 8. (a) Observed ASO season climatology of vertical shear of zonal wind (unit: m s' 1 ) 
between 200 and 850 hPa and mean bias in the (b) FSU, (c) GFDL, (d) GISS, (e) GSFC, and (f) 
GFS models, as well as (g) the MME. The boxes with dash lines denote the MDR. 

Fig. 9. Composites of ASO seasonal mean vertical wind shear anomalies (unit: m s' 1 ) for 
(a),(d) EP El Nino, (b),(e) CP El Nino, and (c),(f) La Nina during 1982-2009 in observations 
(top row) and the MME mean (bottom row). The anomalies circled by light lines are above the 
90% significance level. The boxes with dash lines denote the MDR. 

Fig. 10. Composites of ASO seasonal mean precipitation anomalies (unit: mm day' 1 ) for 
(a),(d) EP El Nino, (b),(e) CP El Nino, and (c),(f) La Nina during 1982-2009 in observations 
(left column) and the MME mean (right column). The anomalies circled by light lines are above 
the 99% significance level. 

Fig. 11. Changes in vertical wind shear (unit: m s' 1 K' 1 ) associated with a westward shift 
of warm SST anomalies from the Nino-3 region (5°S - 5°N, 90° - 150°W) to the Nino-4 region 
(5°S - 5°N, 150°W - 160°E). The boxes with solid lines denote the MDR. 
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598 Table 1. List of five GCMs for the HWG interannual experiments, the number of 

599 ensemble members, model data grid, and references for TC tracking algorithms. 

600 


Model 

Ensemble 

members 

Model data grid 
zonalx meridional 

Tracking algorithm 

FSU 

3 

384 x 192 

LaRow et al. (2008) 

GFDL 

3 

576 x 360 

Zhao et al. (2009) 

NASA GISS 

3 

360 x 180 

Camargo and Zebiak (2002) 

NASA GSFC 

2 

576 x 361 

LaRow et al. (2008) 

NCEP GFS 

5 

360 x 181 

Camargo and Zebiak (2002) 


601 

602 
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603 Table 2. List of TC statistics for observations (OBS), multiple model ensemble (MME) 

604 mean, and individual model ensemble means, including 28-yr (1982-2009) long-term mean 

605 annual number of Atlantic TCs, variance of interannual variability, linear trend (increase of TCs 

606 per decade), anomaly correlation (AC) between observations and model simulated interannual 

607 TC anomalies, and root-mean-square error (RMSE). The variance for each model is the average 

608 of the variance derived from individual ensemble members. 

609 


Model 

Mean 

Variance 

Trend 

AC 

RMSE 

OBS 

11.7 

25.9 

3.7 



MME 

13.1 

17.0 

1.9 

0.86 

3.5 

FSU 

13.5 

9.2 

1.7 

0.62 

4.5 

GFDL 

12.7 

16.4 

2.2 

0.74 

3.6 

GISS 

6.2 

8.8 

1.1 

0.68 

6.7 

GSFC 

10.9 

24.5 

2.6 

0.62 

4.2 

GFS 

22.0 

26.1 

2.1 

0.73 

10.9 


610 

611 
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612 Table 3. Mean annual number of TCs over the entire 28 years, five EP El Nino, five CP 

613 El Nino, and eight La Nina years, respectively, for observations (OBS), MME, and individual 

614 model ensemble means. Values in parentheses are the percentages of the 28-yr climatology. 

615 


Model 

Mean 

EP El Nino 

CP El Nino 

La Nina 

OBS 

11.7 

6.8 (58%) 

10.2 (87%) 

14.6(125%) 

MME 

13.1 

10.4 (80%) 

12.4 (95%) 

14.5 (111%) 

FSU 

13.5 

11.9(88%) 

12.0 (89%) 

15.3 (113%) 

GFDL 

12.7 

9.5 (75%) 

11.0(87%) 

15.5 (122%) 

GISS 

6.2 

4.5 (73%) 

5.7(91%) 

7.3 (117%) 

GSFC 

10.9 

6.9 (63%) 

12.9(118%) 

11.2 (102%) 

GFS 

22.0 

19.3 (88%) 

20.5 (93%) 

23.1 (105%) 


616 

617 


27 


(a) EP El Nino 



120E 150E 180 150W 120W 90W 60W 30W 0 


(c) La Nina 



-1 - 0.75 - 0.5 - 0.25 0.25 0.5 0.75 1 1.25 1.5 


618 

619 Fig. 1. Composites of ASO seasonal mean SST anomalies (unit: K) for (a) EP El Nino 

620 (1982, 1986, 1991, 1997, and 2006), (b) CP El Nino (1987, 1994, 2002, 2004, and 2009), and (c) 

621 La Nina (1983, 1984, 1988, 1995, 1998, 1999, 2005, and 2007) during 1982-2009. The 

622 anomalies circled by light lines are above the 99% significance level estimated by the Monte 

623 Carlo tests. 
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625 

626 Fig. 2. Composites of TC track density (top row) and track density anomaly (middle row) 

627 for (a),(d) EP El Nino, (b),(e) CP El Nino, and (c),(f) La Nina years, and distribution of TC 

628 origins during (g) five EP El Nino, (h) five CP El Nino, and (i) five La Nina years derived from 

629 observations. The anomalies circled by light lines (middle row) are above the 90% significance 

630 level estimated by the Monte Carlo tests. The boxes with dash lines denote the main 

631 development region (MDR; 10°-20°N, 20°-80°W). 
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Fig. 3. Time series of annual number of Atlantic TCs from 1982 to 2009 for observations 
(OBS) and multi-model ensemble (MME) mean (thick lines with open circles), as well as 
individual model ensemble means (thin lines). Grey shading denotes the range of ± one standard 
deviation of the spreads of the five individual model ensemble means around the MME mean. 
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Fig. 4. Climatology of track density for (a) observations, (c)-(g) individual model 
ensemble means, and (h) MME mean, and 28-yr total TC origins for (b) observations, (i)-(m) 
one ensemble member of each model, and (n) MME total from one member of each model. The 
boxes with dash lines denote the MDR. 
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646 

647 Fig. 5. Composites of track density during EP El Nino (left column), CP El Nino (middle 

648 column), and La Nina (right column) for five individual model ensemble mean (top five rows) 

649 and MME mean (bottom row). The boxes with dash lines denote the MDR. 
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Fig. 6. Composites of track density anomaly during EP El Nino (left column), CP El Nino 
(middle column), and La Nina (right column) for five individual model ensemble mean (top five 
rows) and MME mean (bottom row). The anomalies circled by light lines are above the 90% 
significance level estimated. The boxes with dash lines denote the MDR. 
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Fig. 7. Distribution of TC origins during five EP El Nino (left column), five CP El Nino 
(middle column), and five La Nina (right column) years from one ensemble member of each 
model (top five rows) and MME total from one member of each model (bottom row). The boxes 
with dash lines denote the MDR. 
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Fig. 8. (a) Observed ASO season climatology of vertical shear of zonal wind (unit: m s' 1 ) 
between 200 and 850 hPa and mean bias in the (b) FSU, (c) GFDL, (d) GISS, (e) GSFC, and (f) 
GFS models, as well as (g) the MME. The boxes with dash lines denote the MDR. 
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668 

669 Fig. 9. Composites of ASO seasonal mean vertical wind shear anomalies (unit: m s' 1 ) for 

670 (a),(d) EP El Nino, (b),(e) CP El Nino, and (c),(f) La Nina during 1982-2009 in observations 

671 (top row) and the MME mean (bottom row). The anomalies circled by light lines are above the 

672 90% significance level. The boxes with dash lines denote the MDR. 
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674 

675 Fig. 10. Composites of ASO seasonal mean precipitation anomalies (unit: mm day' 1 ) for 

676 (a),(d) EP El Nino, (b),(e) CP El Nino, and (c),(f) La Nina during 1982-2009 in observations 

677 (left column) and the MME mean (right column). The anomalies circled by light lines are above 

678 the 99% significance level. 
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680 

681 Fig. 11. Changes in vertical wind shear (unit: m s' 1 K" 1 ) associated with a westward shift 

682 of warm SST anomalies from the Nino-3 region (5°S - 5°N, 90° - 150°W) to the Nino-4 region 

683 (5°S - 5°N, 150°W - 160°E). The boxes with solid lines denote the MDR. 
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